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NUMERICAL  SIMULATIONS  OF  THE  WAKE  DOWNSTREAM 
OF  A  TWIN-SCREW  DESTROYER  MODEL 


1.  INTRODUCTION 

The  purpose  of  this  report  is  to  present  and  discuss  the  results  of  numerical  simula¬ 
tions  of  high  speed  surface  ship  model  experiments  which  have  recently  been  conducted 
at  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center  (DTXSRDC). 
In  these  experiments  the  turbulent  wake  behind  a  self-propelled  model  of  a  twin-screw 
destroyer  was  documented.  Data  from  two  of  the  experiments  were  forwarded  to  the 
Naval  Research  Laboratory  (NRL)  for  the  purpose  of  determining  the  accuracy  with 
which  such  flow  fields  can  be  numerically  simulated  by  means  of  three-dimensional 
marching  codes.  These  data  were  used  to  generate  initial  conditions  for  calculations 
with  the  NRL  code  TWAKE.  In  the  sections  that  follow  brief  descriptions  of  the  exper¬ 
iments  and  the  procedures  for  initializing  the  computations  are  given.  The  results  of 
the  simulations  are  then  presented  and  discussed.  Comparisons  of  the  calculations  and 
experiments  have  not  been  presented  herein  pending  the  release  of  the  experimental 
results. 

2.  DESCRIPTION  OF  EXPERIMENTS 

Experiments  were  conducted  to  determine  the  three  dimensional  turbulent  velocity 
field  and  the  free  survace  elevations  in  the  wake  region  of  a  high  speed  surface  ship 
model.  A  line  drawing  of  the  model  is  contained  in  Fig.  1  and  Table  1  displays  relevant 
data  on  the  model  dimensions  and  towing  conditions. 

A  rectangular  coordinate  system  fixed  to  the  model  was  chosen  to  define  the  mea¬ 
surement  locations.  The  origin  of  the  coordinate  system  is  at  the  intersection  of  the 
undisturbed  free  surface,  the  ship  aft  perpendicular  (AP),  and  the  vertical  plane  run¬ 
ning  down  the  ship  longitudinal  centerline  (Fig.  2a).  The  “A'”  or  streamwise  axis  is 
directed  aft  (positive  direction)  from  the  origin  and  in  the  model  centerline  plane.  The 
aYn  or  transverse  axis  is  directed  to  starboard  from  the  origin.  The  "Z"  or  vertical 
axis  is  directed  upward  to  complete  the  orthogonal,  right-handed  coordinate  system. 

Mean  velocity  and  turbulence  measurements  were  taken  using  both  laser  floppier 
velocimetry  (LDV)  and  hot-film  annemometry  (HF)  techniques.  Details  of  the  instru¬ 
mentation,  procedures,  and  data  reduction  may  be  found  in  Ref.  [1],  The  measurement 
domain  consisted  of  a  series  of  vertical  (Y  —  Z)  planes  at  various  axial  locations  in  the 
wake.  A  schematic  of  the  typical  experimental  cross-plane  including  the  major  wake 
producing  elements  of  the  model  is  shown  in  Fig.  2b.  Generally  the  sensors  were  placed 
at  a  specified  depth  (Z)  and  measurements  were  taken  point  by  point  at  intervals  of 
2.0  in  in  the  range  —3.0  <  Y  <  33.0.  The  measurement  interval  in  the  Z  direction 
was  also  2.0  in  from  —1.6  in  (nearest  to  the  free  surface)  to  a  maximum  depth  of  be¬ 
tween  —  17.0  in  to  — 2S.0  in  depending  on  the  distance  aft  of  the  model.  At  a  few  axial 
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locations  the  port  half-plane  w as  measured  to  check  the  symmetry  of  the  How  field. 
The  LDV  tests  yielded  data  for  the  three  components  of  mean  velocity,  (£/,Vr, IF),  and 
the  three  turbulent  normal  stresses,  which  can  be  summed  to  give  the  total  turbulence 

kinetic  energy, 


( u'u '  +  v'v'  4-  w'w')  . 


(1J 


The  LDV  data  were  corrected  for  the  wave  orbital  velocities.  The  cross-film  sensors 
were  employed  in  two  orientations  to  yield  independent  realizations  of  the  above  mean 
flow  and  turbulence  quantities  as  well  as  the  two  principal  Reynolds  stresses  ( u'v',u'w'). 

A  matrix  showing  the  scope  of  the  measurements  obtained  is  given  in  Table  2. 
The  numbers  in  the  table  refer  to  the  approximate  number  of  points  measured  for  a 
particular  test  which  is  characterized  by 

1)  the  distance  from  AP  to  the  measurement  plane, 

2)  the  tow  speed, 

3)  the  rotation  direction  of  the  propellers  (or  unpropelled),  and 

4)  the  type  sensor  used  for  the  measurements. 

The  elements  denoted  by  an  asterisk  have  been  forwarded  to  NRL  via  Refs.  [2-3]. 

Two  of  the  data  sets  referenced  in  Table  2  have  been  used  to  develop  initial 
conditions  for  the  numerical  simulations  of  the  evolving  wake  flow  fields.  These  are 
the  4.0  knots  data  at  the  10.0  ft  plane  with  both  outboard  (clockwise)  and  inboard 
(counter-clockwise)  propeller  rotation.  The  initialization  process  for  TWAKE  and  a 
brief  description  of  the  simulation  model  are  given  in  the  next  section. 


3.  INITIALIZATION  OF  TWAKE 


The  NRL  code  TWAKE  is  a  suitably-modified  derivative  of  the  3DPNS  finite- 
element  code  described  in  Ref.  [4].  The  code  solves  the  steady,  three-dimensional,  time- 
averaged  (in  the  turbulence  sense)  parabolic  Navier-Stokes  equations.  The  effects  of 
turbulence  are  described  using  modelled  transport  equations  for  the  turbulence  kinetic 
energy  and  the  isotropic  dissipation  function,  e,  along  with  an  anisotropic  closure  for 
the  turbulent  stresses.  For  later  reference  the  appropriate  equation  set  is  given  below: 


continuity:  __ 
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axial  momentum: 
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vertical  momentum: 


transverse  momentum: 
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In  its  standard  form  the  model  contains  the  eight  constants  {C’i.  C-j.  Ct.  C’t.  - 
C«!,  C( ,}  which  take  on  the  commonly  accepted  values  {0.94.  0.0G7.  0.00.  0.0(1 
1.3.  1.44,  1.92}.  Cooper  in  Ref.  [o]  has  consul  ered  variations  in  these  values  fo 
near  to  a  free  surface.  No  adjustments  were  undertaken  in  the  present  studv. 

The  above  equation  set  describes  an  initial  value  problem  that  requires  di 
tions  of  the  three  components  of  mean  velocity,  the  turbulence  kinetic  eueruv.  a 
dissipation  function  in  order  to  begin  a  computation.  F<»r  the  computations  ;>n 
in  this  report  the  domain  of  Fig.  2b  was  modified  to  include  the  square  iv::mi  - 
bv  0.0  <  1  <  30, 0  m  and  —30.0  <  Z  <  0.0  in.  In  this  domain  triangular  Unite  ei. 


were  constructed  by  connecting  computational  nodal  points  placed  at  1.0  in  intervals 
in  both  the  Y  and  Z  directions.  This  resulted  in  1369  nodes  and  2592  finite  elements, 
a  discretization  approximately  twice  as  dense  as  that  used  in  the  experiments. 

For  each  simulation  (outboard  and  inboard  rotations  at  4.0  knots  beginning  at  the 
X  =  10.0  ft  plane)  the  LDV  data  for  the  mean  velocity  components  and  the  vertical 
and  streamwise  turbulence  normal  stresses  were  linearly  interpolated  onto  the  compu¬ 
tational  grid.  Due  to  the  inherent  difficulty  in  measuring  the  on-axis  (Y)  component 
of  fluctuating  velocity  with  the  LDV  system,  the  data  for  v'v'  was  not  used.  Rather, 
the  turbulence  kinetic  energy  was  formed  by  adding  twice  the  vertical  normal  stress  to 
the  streamwise  normal  stress.  The  dissipation  function  was  assumed  to  be  related  to 
the  turbulence  kinetic  energy  through  a  dissipation  length  scale, 


The  outboard  rotation  HF  tests  at  10.0  ft  provided  data  for  the  principal  components 
of  Reynolds  stress  u'v'  and  u'w'.  Using  Eq.  (14),  the  leading  terms  of  Eqs.  (9)  and 
(10)  were  calculated  from  the  interpolated  experimental  data  for  a  range  of  values  of  Id 
and  compared  to  the  experimental  stresses.  It  was  found  that  a  value  Id  =  0.5S  ft  was 
sufficient  to  maintain  the  dominant  character  of  the  experimental  stress  distributions, 
and  to  preserve  the  characteristic  magnitude  of  these  stresses.  This  value  corresponds 
to  approximately  85%  of  the  propeller  diameter,  23%  of  the  beam,  and  71%  of  the 
draft.  Since  the  inboard  rotation  tests  did  not  include  hot-film  measurements,  the 
same  value  of  Id  was  used  in  the  simulation  of  those  experiments.  Dependent  variables 
at  computational  points  located  outside  of  the  measurement  domain  were  smoothly- 
blended  to  values  consistent  with  the  boundary  conditions,  freestream  conditions  at 
the  bottom  and  right  boundaries  and  symmetry  conditions  at  the  Z  =  0.0  and  Y  =  0.0 
planes.  In  the  following  section  the  results  of  the  simulations  are  discussed. 

4.  RESULTS  OF  THE  SIMULATIONS 


TWAKE  solves  the  governing  equation  system  in  a  non-dimensionalized  form 
wherein  the  independent  variables  are  non-dimensionalized  by  a  characteristic  length 
and  the  independent  variables  are  non-dimensionalized  by  the  characteristic  density 
and  appropriate  powers  of  the  characteristic  velocity.  These  were  taken  as  the  model 
ship  beam,  B  —  2.5  ft,  the  ambient  density,  po  —  1-935  /6/-sec2-/t-4,  and  the  tow  ve¬ 
locity,  Uo  =  6.76  /f-sec-1.  The  origin  for  the  X  axis  was  further  offset  by  the  distance 
from  the  AP  to  the  initial  computational  plane.  A*o  =  10.0  ft  in  both  simulations. 
In  the  following  sub-sections,  for  each  of  the  simulations  in  turn,  full  planar  numer¬ 
ical  solutions  (for  the  experimentally  observed  fluid  properties)  are  presented  at  the 
computational  planes  corresponding  to  the  measurement  locations  shown  in  Table  2. 

Outboard  Propeller  Rotation 

The  first  entries  in  each  of  Figs.  3  and  4  are  contour  plots  of  the  non-dinieusiouni. 
interpolated,  experimental  data  at  the  X  =  10.0  ft  station  for  the  streamwise  mean  ve¬ 
locity  and  the  turbulence  kinetic  energy,  respectively.  The  initial  velocity  distribution 
is  characterized  by  a  drag  wake  (  U /U»  <  1.0)  near  the  free  surface'  and  an  overt lnu.-ted 


region  (U /Uo  >  1.0)  in  the  vicinity  of  the  propeller  disk  (projected  aft).  The  over- 
thrusted  fluid  is  bifurcated  into  distinct  regions,  each  having  local  maximum  velocities 
of  approximately  O.OSUq.  The  maximum  velocity  deficit  in  the  drag  wake,  which  also 
exhibits  several  local  extrema,  is  0.00&'o.  The  initial  kinetic  energy  distribution  con¬ 
tains  local  maxima  in  the  propeller  disk  and  immediately  at  the  free  surface.  The 
maximum  magnitude  corresponds  to  a  characteristic  velocity  of  turbulence  exceeding 
0.06Uo-  The  initial  distribution  of  the  non-dimensional  swirl  velocity, 

va  =  ±(V2+w2y'\ 

Uo 

is  shown  is  Fig.  5a.  There  are  two  distinct  vortical  structures  present  which  are 
approximately  coincident  with  the  two  overthrust  lobes.  Whereas  the  overthrust  is 
about  the  same  in  both  lobes  the  two  structures  are  quite  different.  The  uppermost 
and  strongest  is  evidently  due  to  the  propeller  and  the  swirl  magnitude  approaches 
O.lOUo-  The  lower  counter-rotating  vortex  presumably  was  shed  from  the  hull  or  the 
control  surfaces.  The  fluid  in  the  propeller  swirl  is  considerably  more  turbulent  than 
that  in  the  lower  vortex  (Fig.  4a).  At  the  depth  Z/B  %  —0.4,  the  two  vortices  are 
mutually  acting  to  drive  fluid  toward  the  ship  centerline.  This  action  is  assumed  to  be 
the  cause  of  the  bifurcation  shown  in  Fig.  3a.  Figure  6a  is  the  initial  distribution  of 
the  resultant  of  the  two  principal  Reynolds  stresses, 

Tz  =  -T2-(( U'V ')2  -f  (u'w')2)1/2. 

u  0 

In  this  case  the  displayed  data  are  not  the  actual  (interpolated)  measurements  but 
rather  are  derived  from  Eqs.  (9-10)  and  (14)  using  the  experimental  mean  velocity  and 
turbulence  fields.  The  measured  characteristic  level  of  shear  has  been  preserved  through 
the  choice  of  the  dissipation  length  scale.  The  computed  stress  is  predominantly  radial 
in  the  overthrusted  fluid  and  has  a  dominant  vertical  component  in  the  drag  wake. 

Figures  3b  through  3d  are  contour  plots  of  the  computed  mean  streamwise  velocity 
at  each  of  the  three  downstream  measurement  planes  located  16.0,  22.0,  and  30.0  ft 
aft  of  the  model,  respectively.  The  contours  near  the  free  surface  (drag  region)  show 
that  by  X  =  30.0  ft  the  several  local  extrema  in  the  velocity  have  diffused  to  a 
single  extremum  near  Y/B  =  0.3.  As  the  solution  proceeds  downstream  the  minimum 
velocity,  or  maximum  drag  velocity  at  the  surface,  moves  to  outboard,  existing  finally  at 
Y/B  =  0.4.  The  overthrusted  portion  of  the  velocity  distribution  in  the  subsurface  wake 
essentially  maintains  the  bifurcated  structure  characteristic  of  the  initial  condition.  The 
maximum  velocity  peaks  are  convected  inboard  by  the  propeller  swirl  and  finally  exist 
near  the  model  centerline.  The  computation  shows  that  by  X  =  30.0  ft.  overt  hr  tasted 
fluid  (U /Uo  >  1-0)  has  been  carried  to  the  free  surface  by  the  secondary  flow.  At 
this  downstream  plane  the  maximum  subsurface  overthrust,  is  approximately  0.04t 
while  the  maximum  at  the  surface  is  0.02f'().  The  profiles  of  the  evolving  turbulence 
kinetic  energy  (Figs.  4b-4d)  also  exhibit,  the  convective  action  of  the  swirl  flow.  As  the 
solution  proceeds  downstream  the  maximum  values  of  turbulence  energy  e>:i>t  at  t la- 
free  surface  near  Y/B  =  0.2. 

The  downstream  decay  of  the  swirl  velocity  distribution  is  observed  in  Fig--.  5b 
through  5d.  At  A”  =  30.0  Jt  the  swirl  litis  decayed  to  about  one-third  * > t  the  initial 
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level.  The  lower  inboard  rotating  vortex,  initially  existing  in  fluid  of  considerably  less 
turbulence  than  the  upper,  decays  at  a  relatively  slower  rate  such  that  by  X  =  30.0  ft 
the  two  vortices  are  of  comparable  strength.  The  evolving  shear  stress  distributions 
shown  in  Figs.  6b  through  6d  show  that  by  the  final  measurement  station,  the  stress 
is  predominantly  horizontal  near  the  free  surface.  The  shear,  of  course,  corresponds 
directly  with  the  velocity  gradient  field  apparent  from  the  contour  spacing  of  Figs.  3a 
through  3d. 

Figures  7a  and  7b  are  contour  plots  of  the  streamwise  velocity  and  turbulence 
kinetic  energy,  respectively,  in  the  plane  of  the  mean  free  surface  (Z / B  =  0.0).  In  these 
figures  the  port-side  of  the  wake  ( Y/ B  <  0.0)  has  been  depicted  by  symmetric  reflection 
of  the  computed  starboard-side  solution.  The  frame  of  reference  for  the  velocity  has 
been  transformed  to  a  stationary  one  such  that  negative  contour  values  depict  fluid 
moving  in  the  direction  of  ship  motion  (left  in  Fig.  7a)  and  positive  contour  levels 
denote  overthrusted  fluid  moving  in  the  opposite  direction.  The  contours  show  that 
the  three  velocity  minima  on  each  side  of  the  wake  merge  to  a  single  minimum  by 
{X  —  Xo)/B  as  2.0  ( X/B  =  15.0  ft)  and  that  this  minimum  slowly  moves  away  from 
the  ship  centerline.  The  minimum  velocity  (maximum  drag  velocity)  at  the  end  of  the 
computational  domain  is  about  4%  of  the  tow  velocity.  Overthrusted  fluid  from  the 
propeller  thrust  begins  to  reach  the  surface  at  ( X  —  Xq )/B  =  5.0,  which  corresponds 
to  nearly  18.0  ft  aft  of  the  model.  As  the  wake  turbulence  decays  in  the  downstream 
direction  (Fig.  7b),  the  peak  values  in  the  surface  plane  move  toward  the  ship  centerline. 
The  characteristic  magnitude  of  the  fluctuating  turbulence  velocity  defined  as  A'0-5  is 
approximately  0.0417o  at  the  end  of  the  computation.  This  is  of  the  same  order  as  the 
maximum  mean  velocity  in  the  drag  wake. 

Figures  8a  and  Sb  are  plots  of  the  longitudinal  decay  of  certain  characteristic  values 
of  the  streamwise  velocity  deficit  (or  excess)  and  turbulence  kinetic  energy,  respectively. 
The  curves  labeled  “maximum”  refer  to  the  variation  of  the  maximum  value  of  the  vari¬ 
able  anywhere  in  the  cross-plane  ( X/B  =  const.)  whereas  those  labeled  “maximum  in 
surface  plane”  refer  to  the  variation  of  the  maximum  value  of  the  variable  in  the  sur¬ 
face  plane  (Z/5  =  0.0)  only.  For  the  velocity  the  variation  of  the  minimum  (maximum 
deficit)  is  also  shown.  While  the  curves  are  generally  self-explanatory,  Fig.  Sb  deserves 
some  comment.  This  figure  and  Fig.  4a  show  that  essentially  identical  local  max¬ 
ima  in  the  kinetic  energy  exist  in  both  the  surface  plane  and  the  subsurface  region. 
However,  after  a  developmental  region,  the  turbulence  in  the  surface  plane  decays  at 
a  faster  rate.  Still  further  along  in  the  developing  wake  the  trend  reverses  and  the 
subsurface  turbulence  decays  at  a  faster  rate.  Finally  at  the  end  of  the  computation 
the  maximum  turbulence  is  at  the  free  surface.  This  behavior  can  be  explained  with 
recourse  to  several  of  the  preceeding  figures.  Figures  3a  and  6a  show  that  in  the  initial 
plane  the  larger  values  of  the  Reynolds  stresses  and  the  velocity  gradients  exist  in  the 
subsurface  region.  The  production  term  in  Eq.  (12)  is  therefore  higher  in  this  region. 
The  dissipation  resembles  the  kinetic  energy  (Fig.  4a)  through  Eq.  (14)  and  is  of 
comparable  magnitude  in  both  regions.  Thus  the  balance  is  such  that  the  turbulence 
at  the  surface  will  decay  faster  in  the  upstream  part  of  the  wake.  As  the  wake  develops 
(Figs.  3c  and  6c)  the  production  of  turbulence  becomes  comparable  in  both  regions 
and  thereafter  convection  plays  an  increasingly  important  role  transporting  turbulent 
fluid  to  the  surface.  The  onset  of  surface  renewal  with  fluid  from  below  is  apparent  in 
Fig.  Sa  near  (A'  -  A'o  )/B  =  5.0. 


Inboard  Propeller  Rotation 

The  initial  conditions  for  the  simulation  of  the  inboard  rotation  experiment  are 
shown  as  the  first  entries  in  each  of  Figs.  D  through  12.  Although  the  Fronde  and 
Reynolds  numbers  are  identical  for  the  inboard  and  outboard  rotation  experiments, 
the  measured  fluid  properties  are  significantly  different  at  the  initial  X  =  10.0  ft 
axial  plane.  Figure  9a  shows  that  the  velocity  varies  from  O.STt'o  hi  the  drag  region 
to  l.llC/o  in  the  propeller  thrust  These  extrema  are  significantly  different  than  the 
limiting  values  of  O.OlL^o  and  l.OSt/o  observed  in  the  outboard  rotation  experiment. 
Also  of  note  in  this  figure  is  the  presence  of  a  single  velocity  peak  in  the  thrusted  fluid 
located  relatively  more  to  the  outboard,  rather  than  the  more  inboard  and  bifurcated 
structure  manifest  in  the  first  experiment.  Near  the  free  surface  at  the  model  centerline 
there  is  a  small  region  of  excess  velocity  ( U /Uq  ~  1.02).  Although  this  small  velocity 
excess  may  be  comparable  to  the  experimental  error,  the  qualitative  structure  of  the 
isovels  in  this  area  suggests  a  distinctly  weakened  or  even  non-existent  drag  wake  along 
the  ship  centerline.  The  initial  distribution  of  turbulence  kinetic  energy  (Fig.  10a) 
shows  maximum  levels  50%  greater  than  in  the  outboard  experiment.  In  this  case 
the  fluid  near  the  free  surface  is  considerably  less  turbulent  than  that  in  the  propeller 
swirl.  The  initial  propeller  swirl  (Fig.  11a)  is  considerably  less,  only  about  0.0G3t’o- 
and  located  more  to  outboard  which  is  consistent  with  the  velocity  profile  in  Fig.  9a. 
The  higher  thrust  velocity  and  lower  swirl  suggests  more  efficient  propeller  operation 
in  this  experiment.  The  lower  vortex  apparent  in  the  outboard  data  is  absent.  There 
is  a  significant  outboard  rotating  vortex  near  the  free  surface.  This  vortex,  possibly 
induced  by  the  propeller  swirl,  may  be  responsible  for  transporting  overthrusted  fluid 
to  the  surface  and  producing  the  weakened  drag  wake  noted  in  the  discussion  of  Fig. 
9a.  The  initial  Reynolds  stress  distribution  is  shown  in  Fig.  12a  and  is  characterized 
by  predominantly  radial  shear  in  the  thrusted  fluid.  The  magnitude  is  a  factor  of  two 
higher  than  in  the  first  experiments  due  to  the  higher  turbulence  kinetic  energy  and 
mean  velocity  gradients. 

Table  2  shows  that  data  was  collected  during  the  inboard  rotation  tests  at  two 
additional  downstream  planes,  ,Y  =  16.0  ft  and  .Y  =  30.0  ft.  Although  these  particular 
tests  were  made  while  towing  the  model  at  6.0  knots  rather  than  at  the  speed  of  4.0 
knots  used  in  the  measurements  at  .Y  =  10.0  ft,  it  may  be  possible  to  compare  the 
results  of  the  simulation  to  this  data  in  a  qualitative  manner.  The  second  and  third 
entries  in  each  of  Figures  9  through  12  are  the  computed  results  for  the  selected  fluid 
variables  at  the  two  downstream  stations. 

Figures  9b  and  9c  show  that  as  the  wake  evolves  to  X  =  30.0  ft  the  surface 
drag  velocity  decreases  to  about  0.03bro  while  the  maximum  propeller  overthrust  has 
decreased  to  0.04b*o.  Bv  the  final  station  the  drag  wake  has  been  diffused  and  adverted 
into  the  center  of  the  wake.  The  maximum  level  of  turbulence  kinetic  energy  (Fig. 
10c)  has  decayed  to  approximately  one-third  the  magnitude  that  existed  in  the  initial 
plane.  The  maximum  turbulence  levels  remain  well  below  the  plane  of  rite  mean  free 
surface.  In  contrast  to  the  outboard  simulation,  there  is  no  obvious  effect  on  either 
the  mean  velocity  or  turbulent  energy  due  to  convection  by  the  swirl  flow.  The  swirl, 
shown  in  Figs,  lib  and  11c.  decays  to  20%  of  the  initial  level  by  A*  =  30.0  ft.  The 
computed  Reynolds  stress  distributions  for  the  downstream  planes  are  shown  in  Figs. 
12b  and  12c.  The  dominant  horizontal  component  of  shear  stress  which  was  present  in 
the  outboard  experiments  is  absent  by  X  =  30.0  ft. 
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Figures  13a  and  13b  are  contours  plots  of  the  streamwise  velocity  deficit  and  the 
turbulence  kinetic  energy  in  the  plane  of  the  mean  free  surface.  These  distributions 
strongly  resemble  their  outboard  simulation  counterparts  (Fig.  7)  discussed  in  the 
previous  section.  The  small  velocity  excess  mentioned  in  connection  with  Fig.  9a  has 
vanished  by  (.Y  —  Xo)/B  —  1.0.  Also  worthy  of  note  in  Fig.  13a  is  the  significantly 
smaller  horizontal  strain  ( dU /dY)  than  that  evident  in  Fig.  7a. 

Figures  14a  and  14b  are  plots  of  the  longitudinal  variation  of  characteristic  val¬ 
ues  of  the  velocity  and  turbulence  kinetic  energy  and  illustrate  portions  of  previously 
described  results  in  a  different  form. 

5.  COMPARISON  AND  FINAL  REMARKS 

Figures  15a  and  15b  are  included  to  provide  a  means  of  comparison  between  the  two 
simulations  discussed  above.  These  calculations,  performed  subsequent  to  the  original 
computations,  also  extended  the  solution  considerably  further  downstream  as  shown  in 
the  figures.  It  should  be  noted  that  these  solutions  began  to  encounter  tne  freestream 
boundaries  in  the  transverse  plane  at  approximately  (-Y  —  A'o  )/B  =  50.0  and.  conse¬ 
quently,  the  results  beyond  that  region  do  not  reflect  the  actual  physical  background 
environment.  Until  that  point,  however,  all  variables  appear  to  have  achieved  a  power 
law  decay  with  respect  to  the  streamwise  coordinate  as  shown  in  the  figures.  Keeping 
in  mind  that  these  flow  fields  are  considerably  more  complex  than  any  of  the  “clas¬ 
sical”  flows,  the  indicated  asymtotic  behavior  seems  reasonable.  The  major  contrasts 
between  the  two  simulations  are  that  the  inboard  rotation  velocities  decay  at  a  faster 
rate  than  their  outboard  rotation  counterparts.  Consistent  with  that  trend,  the  inboard 
rotation  turbulence  decays  more  slowly  resulting  in  a  correspondingly  higher  turbulent 
diffusivity  over  the  computational  domain.  The  ultimate  decay  is  influenced  by  the 
initial  dissipation  distribution  as  well  as  the  constants  of  the  model.  In  the  absence  of 
data  for  the  shear  stress,  the  dissipation  for  the  inboard  simulation  was  obtained  by 
requiring  the  same  value  of  the  length  scale  that  was  used  in  the  outboard  simulation. 
The  local  turbulence  time  scale,  proportional  to  I\/t,  is  then  computed  to  be  different 
for  the  two  flows. 

A  final  observation  to  be  made  from  Fig.  15b  regarding  the  inboard  rotation 
simulation  is  that  the  maximum  level  of  turbulence  kinetic  energy  do  indeed  exist  at 
the  free  surface  after  (X  —  Xq)/B  20.0.  a  characteristic  that  was  observed  for  the 
outboard  rotation  computation  (Fig.  Sb). 

In  summary,  the  principal  characteristics  of  the  computed  flow  fields  are  that: 

1)  In  the  outboard  simulation,  the  propeller  thrust  is  bifurcated  into  two  distinct 
regions  which  migrate  toward  the  centerline  of  the  wake  and.  by  symmetry,  merne 
with  those  formed  by  the  port  propeller.  The  upper  of  the  overt hrusted  re-pious 
eventually  is  converted  to  the  free  surface  by  the  secondary  flow  and  this  acr- 
to  displace  the  original  free  surface  drag  wake  to  the  outboard.  This  evolution 
combined  with  the  eventual  arrival  of  maximum  turbulence  kinetic  energy  lev.u 
at  the  surface  results  in  relatively  high  transverse  strain  rate  and  Reynolds  -hear 
stress. 

2)  In  the  inboard  simulation,  the  propeller  thrust  wake  remains  unimode!  and 
tially  maintains  spatial  stationarity.  Initial  and  small  levels  of  velocity  e>;. at 
the  surface  near  the  wake  centerline  quickly  dissipate  as  the  port  and  -ra;  bo.ud 
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drag  wakes  merge  into  a  single  diffuse  velocity  deficit  region.  This  results  in  con¬ 
siderably  smaller  transverse  strain  than  in  the  outboard  rotation  simulation.  The 
significantly  higher  levels  of  turbulence  kinetic  energy  produce  correspondingly 
greater  values  for  the  turbulent  diffusivity  and  therefore  faster  decay  rates  for  the* 

fluid  properties. 

The  computed  differences  in  results  for  the  two  cases  are  significant  and  stem  from 
differences  in  the  initial  plane  data  together  with  the  assumption  of  e<jual  dissipation 
scales  for  the  two  situations.  In  lieu  of  further  guidance  from  detailed  experiments,  the 
latter  assumption  appears  appropriate.  It  is  fascinating  that  the  differences  seen  in  the 
two  measured  data  sets  may  be  attributed  soley  to  the  change  in  direction  of  propeller 
rotation.  Undoubtably,  further  experimentation  will  be  needed  to  fully  explore  these 
changes  in  a  complex  flow  field.  It  is  highly  recommended  that  these  experiments 
include  provisions  for  flow  visualization  of  the  major  flow  features  and  of  the  surface 
flow  patterns  which  clearly  are  sensitive  to  the  initial  wake  structure. 
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Table  1.  Model  bliip  Parameters 


Model 

DTNSRDC  No. 


High  Speed  Surface  Ship 
54 15-1 


Length* 

18.8  ft 

Beam 

29.9  in 

Draft  9.8  in 

Prop  Dia. 

SAG  in 

Block  Coef. 

0.506 

Tow  Speed 

4.0  knots 

Froude  No. 

0.28 

(length  between  forward  and  aft  perpendiculars) 


Table  2.  Points/data  planes  from  LDV  and  HF  anemometrv 
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Outboard,  LDV  96  128_  'ill 

Outboard.  LDV 


■IH 


Feet  from  AP  -1.7  -1.1  -0.4  6.1 


Speed,  knots  4  6  4  6  4  6 


Outboard,  HF 


Inboard,  LDV 


No  Prop,  LDV  117 


*  (forwarded  to  NRL) 

f(denotes  plane  where  symmetry  about  vertical  plane  was  tested) 


175*  I  138  2S0*t  279 1  200 
151 


73* 


0.50  0.75 
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(a)  X  =  10.0  ft 


(b)  X  =  1G.0  ft 

Fig.  3  -  Contours  of  mean  strenimviso  velocity:  outboard 
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Fig.  6  (Cont'd)  —  Reynolds  stress  distribution:  outboard  rotation 
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Fig.  12  (Cont’d)  —  Reynolds  stress  distribution:  inboard  rotation 
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